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NUMERICAL  STUDY  OF  WAVE  PROPAGATION  IN  A  NON-UNIFORM  FLOW 

ALEX  POVITSKY* 

Abstract.  The  propagation  of  acoustic  waves  originating  from  cylindrical  and  spherical  pulses,  in  a 
non-uniform  mean  flow,  and  in  the  presence  of  a  reflecting  wall  is  investigated  by  Hardin  and  Pope  approach 
using  compact  approximation  of  spatial  derivatives.  The  2-D  and  3-D  stagnation  flows  and  a  flow  around 
a  cylinder  are  taken  as  prototypes  of  real  world  flows  with  strong  gradients  of  mean  pressure  and  velocity. 
The  intensity  and  directivity  of  acoustic  wave  patterns  appear  to  be  quite  different  from  the  benchmark 
solutions  obtained  in  a  static  environment  for  the  same  geometry.  The  physical  reasons  for  amplification  and 
weakening  of  sound  are  discussed  in  terms  of  dynamics  of  wave  profile  and  redistribution  of  acoustic  energy 
and  its  potential  and  kinetic  components.  For  an  acoustic  wave  in  the  flow  around  a  cylinder,  the  observed 
mean  acoustic  pressure  is  approximately  doubled  (upstream  pulse  position)  and  halved  (downstream  pulse 
position)  in  comparison  with  the  sound  propagation  in  static  ambient  conditions. 

Key  words,  aeroacoustics,  non-uniform  mean  flow,  higher-order  compact  scheme,  amplification  of 
sound,  stagnation  flow,  flow  around  cylinder 

Subject  classification.  Fluid  Mechanics 

1.  Introduction.  This  study  investigates  numerically  the  influence  of  strong  mean  flow  gradients  on 
the  directivity  and  strength  of  sound  waves  propagating  in  such  a  flow. 

The  goal  of  this  study  is  to  get  computational  insight  into  physical  mechanisms  of  angular  redistribution 
of  acoustic  energy  and  acoustic  pressure.  The  considered  mean  flows  include  potential  flow  around  a  circular 
cylinder,  the  2-D  and  3-D  stagnation  flows,  and  a  90°  corner  flow.  These  idealized  mean  flows  mimic  real 
world  flows  in  areas  of  strong  sound  reflection  and  generation  such  as  a  leading  edge  of  a  wing  or  a  turbine 
blade,  a  protruding  corner  of  a  wall  cavity,  a  wing-fuselage  intersection,  an  impingement  area  of  a  jet,  and 
a  flow  behind  a  bluff  body. 

The  amplification  or  weakening  of  sound  propagating  in  non-uniform  flows  has  attracted  considerable 
attention  of  researchers.  However,  available  computational  studies  are  mainly  restricted  to  one-dimensional 
mean  flows. 

For  one-dimensional  flow  in  a  duct  with  variable  cross-section,  the  acoustic  pressure  appears  to  con¬ 
centrate  towards  the  duct  centerline  or  the  duct  wall  for  upstream  and  downstream  sound  propagation, 
respectively  [11],  [15].  Motivated  by  engine  noise,  the  authors  of  these  studies  use  the  wave  envelope  method 
(WEM)  to  study  propagation  of  sound  in  one-dimensional  and  quasi  one-dimensional  mean  flows  This  method 
is  based  on  consideration  of  the  governing  acoustic  variables  as  sums  of  harmonics.  Additional  assumptions 
are  made  about  rapid  axial  oscillations  and  slowly  varying  amplitude  and  phase  corrections  caused  by  the 
non-uniformity  of  the  flow. 

Atassi  and  Grzedzinsky  [2]  considered  propagation  of  unsteady  disturbances  in  flows  around  bodies 
with  a  stagnation  point.  For  incompressible  mean  flows,  the  aeroacoustic  problem  is  formulated  in  terms 
of  an  integral  equation  of  the  Fredholm  type.  The  considered  streamlined  body  and  the  corresponding 
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potential  mean  flow  result  from  superposition  of  a  uniform  flow  and  a  source.  Yet,  additional  assumptions 
and  simplifications  (some  of  them  depend  on  the  type  of  potential  background  flow)  are  needed  to  perform 
integration.  Authors  [2]  present  the  disturbance  pressure  along  the  stagnation  streamline  that  shows  increase 
of  the  disturbance  pressure  for  an  upstream  harmonic  disturbance. 

Hardin  and  Pope  approach  [7]  implies  consideration  of  a  flowfield  as  a  sum  of  the  mean  velocity  and 
a  disturbance  field  and  solution  of  disturbance  variables  in  time  and  space  domain  by  explicit  integration 
in  time  and  use  of  appropriate  spatial  discretization.  The  computational  methodology  used  in  the  current 
study  is  based  on  above  approach. 

In  order  to  use  an  analytical  solution  for  mean  flow,  the  flow  in  the  current  study  is  assumed  to  be 
potential,  i.e.,  inviscid  and  incompressible.  The  speed  of  sound  (relative  to  the  local  mean  flow  velocity)  is 
assumed  to  be  equal  to  unity.  These  assumptions  are  fair  unless  the  Mach  number  exceeds  0.5. 

In  this  paper,  the  relative  orientation  of  the  direction  of  wave  propagation,  mean  pressure  gradient, 
and  the  mean  flow  vector  are  studied  to  show  the  amplification  or  weakening  of  the  acoustic  pressure. 
The  consideration  starts  with  two-dimensional  stagnation  mean  flows,  then  the  approach  is  expanded  to 
three-dimensional  stagnation  mean  flows,  and,  finally,  applied  to  the  inviscid  flow  around  circular  cylinder. 

In  all  considered  cases,  the  obtained  acoustic  pressure  is  compared  with  that  for  the  wave  propagation 
under  static  ambient  conditions  and  for  the  same  geometry  of  surrounding  rigid  boundaries. 

Variation  of  directivity  of  acoustic  pressure  implies  changes  in  directivity  of  acoustic  energy  in  comparison 
with  the  sound  propagation  under  static  ambient  conditions.  This  paper  studies  angular  redistribution  of 
acoustic  energy  and  its  potential  and  kinetic  components  in  non-uniform  mean  flows.  The  acoustic  energy 
is  a  sum  of  potential  and  kinetic  components  where  the  potential  energy  is  closely  related  to  the  averaged- 
in-time  root  mean  square  of  acoustic  pressure  which  is  usually  measured  by  experimentalists  as  an  indicator 
of  noise  intensity. 

Our  recent  study  [12]  shows  that  the  generation  of  sound  waves  originated  from  entropy  sources  is  caused 
by  baroclinic  generation  of  disturbance  vorticity  as  a  result  of  interaction  of  the  mean  flow  pressure  gradient 
and  the  wave  density  gradient.  In  turn,  this  generated  vorticity  causes  the  generation  (or  sink)  of  acoustic 
energy  by  the  source  term  in  the  acoustic  energy  equation.  Order  of  magnitude  analysis  suggests  that  the 
source  term  can  not  be  neglected  in  the  acoustic  energy  equation.  This  study  presents  vorticity  and  energy 
source  terms  obtained  from  numerical  simulation.  For  the  propagation  of  acoustic  waves  originated  from 
an  acoustic  pulse  (as  opposed  to  the  sound  generation  from  entropy  or  vorticity  pulses  [12]),  The  major 
effect  of  the  non-uniformity  of  the  mean  flow  on  the  wave  propagation  (as  opposed  to  the  wave  generation 
from  entropy  or  vorticity  pulses  [12])  is  the  angular  redistribution  of  acoustic  energy  and  kinetic-to-potential 
transform  whereas  the  newly  generated  acoustic  energy  appears  to  be  minor. 

The  article  is  built  as  follows.  In  Section  2,  the  wave  propagation  originated  from  a  cylindrical  2-D 
source  in  2-D  stagnation-type  flows  is  studied  and  numerical  results  about  amplification  of  acoustic  pres¬ 
sure,  alternation  of  wave  pattern,  and  acoustic  energy  redistribution  are  presented.  In  Section  3,  the  wave 
propagation  from  a  spherical  source  is  considered  in  3-D  stagnation  flows.  The  symmetric  3-D  stagnation 
flow  with  a  single  stagnation  point  and  the  3-D  stagnation  flow  with  stagnation  line  flow  are  taken  as  typical 
representatives  of  three-dimensional  stagnation  mean  flows.  In  Section  4,  the  inviscid  flow  around  a  circular 
cylinder  is  considered  for  various  subsonic  Mach  numbers  of  the  mean  flow.  Amplification  of  sound  (for 
upstream  wave  propagation)  and  weakening  of  sound  (for  downstream  wave  propagation)  is  studied  in  this 
Section. 
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2.  Sound  propagation  in  a  2-D  stagnation- type  mean  flow.  The  common  type  of  mean  flow  is 
that  the  flux  across  a  given  surface  is  equal  to  zero,  either  because  the  surface  is  a  symmetry  surface  or 
because  the  surface  is  the  boundary  of  a  rigid  body  [3].  If  two  straight  zero-flux  boundaries  intersecting  at 
an  angle  7r/n  the  stream-function  field  for  incompressible  flow  is  given  by 

=  Arnsin(nO ),  (2-1) 

where  (r,  6)  are  polar  coordinates.  Note  that  ifi  =  0,  if  0  =  0  or  9  =  7 r/n.  For  example,  the  2-D  stagnation 
flow 


U  =  x,V  =  -y 


(2.2) 


is  a  union  of  two  flows  with  n  —  2  where  the  dividing  streamline  is  the  symmetry  axis  of  the  flow. 

Introducing  a  disturbance,  instantaneous  velocities  and  density  are  considered  as  sums  of  the  known 
steady  incompressible  mean  flow  (see  above)  and  the  unsteady  compressible  isentropic  disturbance 

u  =  U  +  uf,  v  =  V  +  vf,p  =  P  +  p/,  p  =  1  +  pt<  (2*3) 


The  normalized  speed  of  sound  is  assumed  to  be  equal  to  unity  and,  therefore,  pressure  and  density  distur¬ 
bances  are  equal  for  acoustic  pulse.  Substituting  above  sums  to  the  Euler  equations,  the  dynamic  equations 
for  unsteady  (disturbance)  components  of  mass  and  momentum  fluxes  are  obtained 


Qt  +  fx+9y  — 


(2.4) 


where 


/  p'  \ 


Q  = 


I  (1  4-  Pf)(U  +  uf)  I 

V  (1  +  pf)(V  +  vi)  ) 


f= 


(1  +  pt)(U  4-  uf) 

(l  +  pf)(U  +  uf)2  +  (P+pt) 
(1  4-  p*){U  +  uf)(V  +  vt) 


(1  +  p,)(V  +  vi)  \ 

(1  +  pf)(U  +  uf)(V  +  vr) 

(1  +p/)(F  +  7;/)2  +  (P  +  p/)  ) 

(2.5) 


The  equations  are  linearized  and  solved  with  respect  to  unknown  disturbance  variables  in  space-time 
domain 


du  _  y 

dt  dx  dx  Ox U  dy 

dv _  dp  y  dv  jjdv  _  dV 

dt  dy  dy  dx  dx 

dp  _  du  dv  ydp  _  y  dp 

dt  ~~  dx  dy  dx  dy' 


dU  ,rdU  TTdU 

_v-V—p-u—p, 

av  TTdv  .. dv 


(2.6) 


The  spatial  derivatives  of  the  mean  flow  variables  for  arbitrary  n  are  presented  in  Appendix  A.  The  non-linear 
set  of  equations  for  the  2-D  stagnation  mean  flow  is  presented  in  our  study  [12]. 


2.1.  Numerical  algorithm  and  the  code  validation.  The  solution  is  advanced  in  time  in  five  sub¬ 
stages  per  time  step  using  a  low-storage  explicit-in-time  fourth-order  Runge-Kutta  (RK)  scheme  scheme 
proposed  by  Williamson  [16]  and  implemented  by  Wilson  et  al.  [17]. 

The  spatial  derivatives  of  unknown  disturbance  variables  are  approximated  using  compact  finite  differ¬ 
ence  schemes  [10]  applied  to  natural  variables  uf ,  vt  and  pt 


+  u'i  +  all'i+1  =  ^^(ui+ 1  -  Ui-i)  +  TK^(Ui+2  -  Ui~2^ 


2Ax 


(2.7) 


3 


where  U[ ,  and  [7^+1  are  unknown  derivatives  at  grid  nodes  i  —  l,i,  and  t  + 1,  Ax  is  the  grid  step  and 

£/j-2j  denote  values  of  variable  U  at  grid  nodes  i  —  2,  ...,i  +  2.  The  values  of  U  are  taken  from  the 

previous  stage  of  the  RK  scheme.  Here  the  classical  fourth-order  Pade  compact  scheme  (a  =  0.25,  a  =  1.5 
and  6  =  0)  which  has  a  tridiagonal  matrix  for  the  right  and  left  sides  of  (2.7)  is  used. 

At  all  boundaries  but  a  rigid  plate  characteristic  inflow  or  outflow  boundary  conditions  are  applied 
to  disturbance  variables.  At  a  rigid  surface,  reflection  boundary  conditions  are  used  for  pressure  distur¬ 
bance  and  for  the  parallel  to  the  surface  component  of  velocity  disturbance  whereas  the  normal  to  the  rigid 
surface  component  of  velocity  is  taken  equal  to  zero.  Discretization  of  spatial  derivatives  in  the  direction 
perpendicular  to  a  boundary  are  computed  by  one-sided  finite  differences  at  all  boundaries  [4]. 

The  system  (2.7)  with  near-boundary  one-sided  differences  is  solved  by  version  of  Gaussian  elimination 
for  narrow  banded  systems  known  as  the  Thomas  algorithm  [9]  applied  to  each  line  of  numerical  grid  in  all 
spatial  directions. 

To  validate  the  code  against  the  known  analytical  solution  [8],  the  static  ambient  conditions  (i.e. ,  the 
mean  flow  velocity  is  equal  to  zero)  are  taken.  In  this  case,  the  right-hand  side  of  system  (2.6)  is  simplified 
as  follows 


du 

dp 

di  = 

dx 

dv 

dp 

dt 

dy 

dp 

du 

dv 

dt 

dx 

dy 

The  initial  2-D  cylindrical  acoustic  pulse  is  given  by 


p  =  p  =  e  exp 


d2  (g  ~  ^c)2  +  {V~  3/c)2 
a 


(2.9) 


where  e  =  0.01, a  —  Jn(2)/9,d  =  60  is  the  normalization  coefficient,  (xciyc)  is  the  pulse  center  coordinates. 

To  check  the  symmetry  of  solution,  the  computational  domain  is  taken  as  a  union  of  [—1, 1]  x  [—1,0] 
and  [-1, 1]  x  [0, 1]  subdomains  separated  by  a  rigid  flat  plate  located  at  y  =  0.  The  couple  of  initial  pulses 
(2.9)  is  located  at  (0,  ±0.25).  Isolines  of  acoustic  pressure  for  this  problem  at  160  x  160  grid  per  subdomain 
are  presented  in  Fig.  1  at  t  =  0.5.  The  isobar  values  are  uniformly  distributed  between  the  maximum  and 
minimum  of  acoustic  pressure.  Although  the  Thomas  algorithm  is,  strictly  speaking,  non-symmetric  (as  any 
Gaussian  elimination),  the  acoustic  pressure  field  in  Fig.  la  is  perfectly  symmetric. 

The  results  of  grid  refinement  study  are  presented  in  Fig.  lb  for  [—1,1]  x  [0, 1]  computational  domain 
covered  with  80  x  80, 120  x  120,  and  160  x  160  numerical  grids.  The  visible  wiggles  exist  only  for  the  coarse 
grid  and  the  solution  approaches  the  analytical  solution  including  most  critical  points  of  the  maximum  and 
minimum  of  the  reflected  wave.  For  the  transmitted  wave,  the  numerical  solution  practically  coincide  with 
the  analytical  solution  [8]  for  finer  grids.  Therefore,  the  grid  160  x  160  has  been  adopted  in  this  study. 

2.2.  Distribution  of  acoustic  pressure.  To  study  the  acoustic  pressure  distribution  in  a  non-uniform 
flow,  the  2-D  stagnation  inviscid  mean  flow  is  taken  as  an  example  (Case  2- A).  The  velocity  field  of  the 
stagnation  mean  flow  and  isobars  of  the  initial  acoustic  pulse,  centered  at  (0,0.25),  are  presented  in  Fig.  2a 
(left).  The  speed  of  sound  is  equal  to  unity  and  is  presented  in  Fig.  2a  for  scale.  As  far  as  the  stagnation 
flowfield  satisfies  Eq.  (2.2),  the  Mach  number  at  a  point  is  equal  to  the  local  distance  from  the  origin,  i.e., 
the  Mach  number  at  the  location  of  the  center  of  initial  pulse  is  equal  to  0.25. 
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The  propagation  of  sound  is  described  by  the  system  (2.6)  where  the  mean  flow  is  given  by  (2.1)  with 
n  =  2.  Grid  refinement  study  for  Case  2- A  is  presented  in  [12]  and  shows  that  the  results  on  the  160  x  160 
and  200  x  200  grids  are  quite  closed  to  each  other  that  confirms  the  use  of  the  former  grid  for  simulations. 

The  isobars  of  acoustic  pressure  are  presented  in  Fig.  2a  (right)  at  t  —  0.5.  Approximate  centerline 
coordinates  of  the  maximums  of  acoustic  pressure  for  the  transmitted  and  reflected  waves  are  0.5  and  0.25, 
respectively.  Note  that  the  Mach  numbers  at  these  points  are  equal  to  their  y-  coordinate  values  (see  above). 
As  expected,  the  acoustic  isobars  obtained  in  the  static  ambient  conditions,  denoted  here  as  Case  2-B,  (Fig. 
la)  are  different  from  those  in  the  Case  2- A  at  the  same  time  moment.  In  the  Case  2- A,  the  y— coordinate 
of  the  wave  front  at  the  neighborhood  of  the  centerline  is  smaller  and  the  x— coordinate  of  the  wave  front 
near  the  wall  is  larger  than  those  in  the  Case  2-B.  The  explanation  is  straightforward  and  is  based  on  the 
fact  that  the  speed  of  sound  in  a  steady  frame  is  equal  to  c  +  £/,  where  c  is  the  local  speed  of  sound  and  U 
is  the  local  velocity  of  the  flowfield. 

The  density  of  isobars  is  larger  at  the  neighborhood  of  the  centerline  and  smaller  near  the  wall  in  the 
Case  2-A  in  comparison  with  those  in  the  Case  2-B.  In  Fig.  3a,  the  acoustic  pressure  profiles  at  the  centerline 
are  presented  for  both  cases  at  t  =  0.5.  The  maximum  absolute  value  of  acoustic  pressure  is  approximately 
80%  larger  for  the  transmitted  wave  and  two  times  larger  for  the  reflected  wave  than  those  in  the  Case  2-B. 
Recall,  that  the  speed  of  wave  propagation  is  c  +  U.  While  the  wave  propagates  upstream  the  decelerating 
subsonic  flow,  its  front  moves  slower  than  its  back.  The  situation  is  opposite  for  the  accelerating  flow  (which 
is  the  case  near  the  wall).  Thus,  the  wave  front  becomes  narrower  at  the  neighborhood  of  the  centerline  and 
wider  while  the  stagnation  flow  spread  along  the  wall  in  comparison  to  those  in  static  ambient  conditions. 

Accelerating  mean  flow  (backward  2-D  stagnation  flow)  is  considered  to  study  the  influence  of  the 
mean  flow  direction  on  the  wave  patterns  (Fig.  2b).  The  initial  position  of  acoustic  pulse,  the  mean  flow 
streamlines  and  the  mean  pressure  distribution  remain  the  same  as  for  the  Case  2-A.  (compare  Fig.  2b  with 
Fig.  2a).  However,  the  direction  of  the  maximum  acoustic  pressure  has  been  considerably  changed  from 
being  perpendicular  to  the  wall  to  being  parallel  to  it.  Still,  the  maximum  acoustic  pressure  corresponds  to 
the  direction  where  the  sound  propagates  upstream  of  the  mean  flow. 

To  study  the  influence  of  the  background  flow  while  the  acoustic  pulse  is  located  off  the  centerline  and 
two  reflecting  surfaces  are  present,  the  90°—  corner  geometry  is  considered  (Fig.  7).  The  rigid  plains  coincide 
with  the  Cartesian  axis.  The  potential  flow  in  this  geometry  is  described  by  a  single  stagnation  flow  with 
n  =  2  flowing  to  the  right.  Initial  position  of  acoustic  pulse  is  taken  (0.25,0.25),  where  the  local  Mach 
number  of  the  mean  flow  is  approximately  equal  to  0.35. 

Computational  results  in  terms  of  acoustic  pressure  are  compared  with  those  obtained  in  the  Case  2-B. 
In  the  latter  case  the  maximum  acoustic  pressure  is  in  the  bisector  direction  and  the  cylindrical  shape  of 
wave  pattern  is  symmetric  with  respect  to  the  bisector  (Fig.  7b).  In  presence  of  the  mean  flow,  the  waves 
flatten,  especially  in  the  neighborhood  of  the  y— axis,  where  the  sound  propagates  upstream  of  the  mean 
flow  (Fig.  7a).  The  maximum  acoustic  pressure  increases  about  40%  in  comparison  with  that  in  the  static 
ambient  conditions. 

2.3.  Distribution  of  acoustic  energy.  To  get  physical  insight  into  amplification  of  acoustic  pres¬ 
sure,  acoustic  energy  is  considered  here.  Transport  equation  of  acoustic  energy  density  for  irrotational  and 
incompressible  mean  flow  and  adiabatic  acoustic  disturbance  is  given  by  [6]: 

^  +  VI  =  u/-(«y  X  <*>'),  (2.10) 

at 
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where 


(2.11) 


E  =  pf2/2  +  u/2/ 2  +  pm/  •  U 

is  the  acoustic  energy  density  and 

I  =  (p,  +  Uf  •  €/)(«/  +  ptU)  (2.12) 

is  the  acoustic  energy  flux. 

The  acoustic  energy  is  a  sum  of  potential  energy  Pot  =  pf2  j 2,  that  mimics  the  strength  of  acoustic 
pressure  in  compressible  media,  and  kinetic  energy  K  =  u/2/2  4*  pm/  •  U. 

The  centerline  distribution  of  E  is  presented  in  Fig.  3b,  where  the  curve  1  corresponds  to  the  Case 
2-A  and  the  curve  2  corresponds  to  the  Case  2-B.  For  the  reflected  wave,  the  local  maximum  of  E  in  the 
Case  2-A  is  about  three  times  greater  than  that  in  the  Case  2-B.  Its  location  is  the  same  as  that  for  the 
maximum  acoustic  pressure  (compare  Fig.  3b  to  Fig.  3a)  and  the  ratio  of  acoustic  energies  is  similar  to  the 
ratio  of  squares  of  acoustic  pressures  (four  times).  It  is  not  the  case  for  the  transmitted  wave,  where  the 
maximum  values  of  E  are  almost  the  same  in  cases  2-A  and  2-B  in  spite  of  the  obvious  difference  in  the 
acoustic  pressure  values. 

To  show  the  total  acoustic  energy  along  the  centerline,  the  integrals 


Ie  = 

f  E(y)dy 

(2.13) 

Jo 

Ip  = 

[  Pot(y)dy 
'0 

(2.14) 

are  plotted  as  functions  of  y  at  the  time  moment  t  =  0.5  in  Fig.  3c.  Obviously,  the  integrals  reach  their 
maximum  when  the  coordinate  y  exceeds  the  front  coordinate  of  the  transmitted  wave.  The  values  of  these 
maximums  for  integral  of  acoustic  energy  differ  only  by  15%  whereas  the  difference  in  integrals  of  potential 
energy  is  about  two  times.  Thus,  the  amplification  of  sound  along  the  centerline  occurs  mainly  by  increase 
of  the  potential  part  in  total  acoustic  energy. 

On  the  contrary,  for  the  bisector  direction  (45°)  and  for  a  near-wall  direction  (75°)  the  part  of  potential 
energy  remains  approximately  the  same  for  the  stagnation  mean  flow  and  for  the  static  ambient  conditions. 
Therefore,  the  weakening  of  sound  near  the  wall  occurs  mainly  due  to  the  angular  redistribution  of  acoustic 
energy.  The  direction-dependent  behavior  of  acoustic  energy  is  explained  by  the  relative  value  of  the  last 
term  in  (2.11).  If  sound  propagates  upstream  of  the  mean  flow,  the  vectors  in  the  inner  product  uf  •  U  are 
collinear,  therefore,  the  variation  of  U  in  a  non-uniform  mean  flow  leads  to  more  prominent  reduction  of  the 
kinetic  energy  in  comparison  with  other  directions.  The  kinetic  energy  transforms  into  potential  energy  and 
leads  to  increase  of  potential  part  in  the  acoustic  energy. 

To  study  the  generation  of  acoustic  energy  by  means  of  interaction  of  the  wave  with  the  non-uniform 
mean  flow,  we  present  here  the  simplified  form  of  the  vorticity  transport  equation  [12]: 

—  =  Vp/  x  VP.  (2-15) 

Dt 

Note  that  the  presence  of  the  pressure  gradient  in  the  mean  flow  (i.e.,  non-uniformity  of  the  flow)  leads 
to  non-zero  right-hand  side  of  above  equation.  Otherwise,  the  vorticity  remains  equal  to  zero  as  in  an 
irrotational  mean  flow.  In  turn,  the  non-zero  ujt  leads  to  the  non-zero  source  term  in  the  energy  equation 
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(2.10).  The  vorticity  is  generated  where  the  acoustic  wave  passes,  i.e.  Vp/  ^  0.  As  fax  as  uf  is  a  multiplier 
in  the  right-hand  side  of  Eq.  (2.10),  the  exchange  of  energy  between  the  steady  mean  flow  and  the  unsteady 
disturbance  occurs  only  in  presence  of  the  wave.  The  generated  vorticity  flows  with  the  local  speed  of 
background  flow  and  eventually  the  vorticity  is  left  behind  the  wave,  that  moves  with  the  sonic  speed.  The 
instantaneous  pattern  of  the  vorticity  generation  source  and  the  acoustic  energy  source  at  t  =  0.5  are  shown 
in  Fig.  5.  As  expected,  areas  of  maximum  acoustic  energy  generation  coincide  with  instantaneous  positions 
of  the  reflected  and  transmitted  waves. 

Approximately,  the  maximum  value  of  the  energy  source  term  is  located  at  the  bisector  direction  (45°). 
In  Fig.  6,  the  acoustic  energy  and  its  source  are  shown  as  functions  of  the  distance  from  the  origin  along  the 
bisector  line  at  t  =  0.5.  The  maximum  value  of  the  source  is  approximately  one  order  of  magnitude  smaller 
than  that  for  the  acoustic  energy.  Therefore,  the  local  gain  of  acoustic  energy  is  small  in  comparison  with 
the  angular  redistribution  of  acoustic  energy  and  kinetic-to-potential  transforms  considered  above. 

3.  Acoustic  pulse  in  the  3-D  stagnation  flow.  In  this  section  we  consider  a  spherical  pulse  prop¬ 
agating  in  a  3-D  stagnation  flow.  The  3-D  stagnation  flow  with  a  single  stagnation  point  at  the  origin  is 
denoted  as  Case  3-A  and  given  by 

W  =  -z,  U  =  0.5®,  V  =  0.5 y.  (3-1) 

The  stagnation  flow  with  the  stagnation  line  (Case  3-B)  is 

W  =  -z,  U  =  x,  V  =  0.  (3.2) 


The  developed  code  is  readily  available  for  3-D  parallel  computations  as  well  [13]. 

In  Case  3-A,  the  propagation  of  disturbance  is  described  by  the  following  system  of  equations 

du 


§i  =  _|P+0.5(. 

dt  dx 


du 

-x- - u  - 

ox 


y 


dy 


dv  dp  _  dv  dv  y.  dv 

dw  dp  .  dw  dw  dw 

=  +  0.5  (-y-^  -  *  )  +*> ~  zP+z-9-z' 


dt 

dp 

dt 


'dy 


'  dx 


op  _  du  dv  dw  h(_xdp  _  v<>Py  ,  zdp 
—  -5T-  —  +a{H  ydy,  +  Zdz ' 


dx  dy  dz  K  dx 

In  Case  3-B,  the  governing  system  of  equations  is  given  by 


du 

dp 

du 

!h  ~ 

dx 

~xdi 

—  u  —  : 

dv 

dp 

dv 

dv 

dt  ~ 

dy 

~xdi 

+  Z!te 

dw 

dw 

i  _ 

~dt  ~ 

_ 

dz 

X  dx 

n  W 

dp  _ 

du 

dv 

dw 

dt 

dx 

dy 

‘  !te  ~ 

du 

xp  +  zTz 


dw 

p  +  zTz' 
<.®L±?®P- 

dx  dz ’ 


(3.3) 


(3.4) 


The  initial  3-D  spherical  acoustic  pulse  is  given  by  an  expression  similar  to  (2.9). 

To  validate  the  code,  the  computations  are  performed  on  a  set  of  numerical  grids  where  the  wave  spreads 
under  static  ambient  conditions  and  in  presence  of  the  rigid  reflecting  plate  (Case  3-C).  The  analytic  solution 
for  acoustic  pressure  in  case  of  the  infinite  domain  is  given  by 

p  -  —  {(r  -  t)  exp[— a(r  -  t)2}  +  (r  +  t)  exp[-fl(r  +  i)2]},  (3.5) 
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where  r  is  the  location  of  the  pulse  center  at  t  =  0.  Solution  for  the  semi-infinite  domain  is  obtained  by  use 
of  the  image  pulse  located  at  —r  : 

p  -p{r)  +p(-r ).  (3.6) 

Results  of  computations  on  803,  1203,  and  1603  grids  are  compared  to  the  above  analytic  solution. 
Profiles  of  the  centerline  acoustic  pressure  on  the  above  set  of  grids  at  t  =  0.5  are  presented  in  Fig.  8.  The 
observations  from  the  grid  refinement  study  are  similar  to  those  obtained  for  the  2-D  case  (Section  2).  The 
numerical  solution  on  the  1603  numerical  grid  practically  coincides  with  the  analytic  solution  and  this  grid 
is  used  for  numerical  simulations. 

To  compare  acoustic  fields  in  Cases  3-A,  3-B,  and  3-C,  the  isolines  of  acoustic  pressure  at  t  =  0.5  are 
shown  in  Fig.  9.  In  Cases  3-A  and  3-C  (Fig.  9a,  b),  the  acoustic  field  is  uniform  with  respect  to  the  direction 
in  the  x  —  y  plane.  The  presence  of  the  stagnation  flow  leads  to  amplification  of  sound  at  the  neighborhood 
of  the  centerline.  In  Case  3-B,  the  wave  pattern  is  direction-dependent:  the  maximum  elongation  of  the 
wave  pattern  along  the  rigid  plane  occurs  in  the  section  x  —  z,  whereas  the  minimum  elongation  is  in  the 
section  y  —  z.  In  the  latter  case,  the  elongation  is  closed  to  that  for  the  Case  3-C  (compare  Fig.  9a  to  Fig. 
9d).  Yet,  the  centerline  position  of  the  pulse  is  the  same  in  Cases  3-A  and  3-B  (see  Fig.  9b, c,  and  d). 

The  centerline  acoustic  pressure  for  Cases  3-A  and  3-B  is  presented  together  with  its  counterpart  for  the 
Case  3-C  at  t-  0.5  in  Fig.  10.  In  spite  of  the  different  wave  pattern  in  Cases  3-A  and  3-B,  the  centerline 
profiles  coincide  for  these  cases.  The  amplification  of  acoustic  pressure  along  the  centerline  is  controlled 
by  redistribution  of  the  potential  and  kinetic  parts  of  acoustic  energy  in  this  direction  as  opposed  to  the 
angular  redistribution  of  acoustic  energy  (see  the  previous  section).  As  far  as  the  centerline  profiles  of  the 
mean  velocity  and  pressure  are  the  same  in  Cases  3-A  and  3-B,  the  profiles  of  acoustic  pressure  appear  to 
be  the  same  along  the  centerline.  Note  that  this  is  quite  a  special  situation  where  both  the  initial  conditions 
of  pulse  and  the  background  flow  velocity  are  the  same  along  the  streamline  for  two  different  flows. 

The  maximum  acoustic  pressure  for  the  transmitted  wave  is  two  times  larger  in  Cases  3-A  and  3-B  than 
that  in  the  Case  3-C.  The  maximum  acoustic  pressure  for  the  reflected  wave  becomes  about  three  times 
larger  in  presence  of  stagnation  flow  than  that  for  the  static  ambient  conditions. 

The  amplification  of  sound  in  the  3-D  case  is  larger  than  that  in  the  2-D  case.  Recall  that  in  both  cases 
the  background  flow  velocity  and  pressure  are  the  same  along  the  centerline.  To  get  qualitative  explanation 
of  such  a  difference,  the  analytic  solutions  for  acoustic  pressure  at  the  centerline  (static  ambient  conditions) 
is  presented  for  the  2-D  and  3-D  cases  in  Fig.  11a.  The  spatial  derivative  of  acoustic  pressure  along  the 
centerline  is  shown  in  Fig.  lib.  By  inspection  of  Fig.  11a  and  by  comparison  of  derivatives  in  Fig.  lib,  one 
can  conclude  that  the  back  and  front  wave  fronts  are  sharper  in  the  3-D  case  than  those  in  the  2-D  case. 
The  front  of  the  wave  propagating  upstream  move  slower  than  the  back  of  it,  that  leads  to  the  amplification 
of  sound  (see  the  previous  section).  The  sharper  the  wave  front  is,  the  more  prominent  is  the  phenomenon. 
Therefore,  the  3-D  pulse  in  a  3-D  stagnation  flow  amplifies  larger  than  that  in  the  2-D  case. 

4.  Aero  acoustics  of  the  flow  over  a  circular  cylinder.  The  propagation  of  sound  waves  from 
a  single  acoustic  pulse  in  the  2-D,  inviscid  irrotational,  incompressible  mean  flow  over  an  infinite  circular 
cylinder  is  considered  here.  The  computational  domain  is  shown  in  Fig.  12.  The  circle  of  unit  diameter 
(Rcyi  —  0.5)  is  the  inner  boundary  and  a  circle  of  3.5  diameter  is  the  outer  boundary.  At  the  outer  boundary 
the  mean  flow  is  practically  uniform.  The  streamlines  of  the  mean  flow  and  the  isobars  of  mean  pressure  are 
computed  analytically  for  potential  flow  [1]  and  presented  in  Fig.  12.  The  pulse  is  located  at  the  centerline, 
in  some  distance  from  the  rigid  cylinder.  While  the  mean  flow  is  directed  from  left  to  right,  the  acoustic 
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waves,  originated  from  the  pulse,  propagate  upstream  of  the  mean  flow.  Otherwise,  the  acoustic  waves 
propagate  downstream  of  the  mean  flow.  The  mean  flow  in  front  of  the  cylinder  can  be  described  by  the 
incompressible  potential  flow  model  fairly  well  up  to  M  =  0.5.  The  after-body  flow  can  be  described  by  this 
model  when  M  <  0.4.  It  has  been  shown  that  the  flow  streamlines  around  cylinder  are  quite  similar  to  those 
for  incompressible  flow  if  M  <  0.4.  [14].  For  higher  Mach  number,  strong  recirculating  zones  appear  at  the 
rear  part  of  the  cylinder  and  the  shock  waves  take  place  at  the  mid-section  of  the  cylinder. 

In  polar  (r,  0)  coordinates,  the  equations  of  propagation  of  small  disturbances  in  a  known  potential  flow 
are  given  by 


du 

dp 

u ^ 

dU  .. 
— —  u  —  U 

~dt  ~ 

dr 

dr 

dr 

dv 

di  ~ 

1  dp 
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dr 1 


do' 


d$ 


d0‘ 


(4.1) 


where  U ,  V,  it,  and  v  are  radial  and  angular  components  of  the  mean  flow  velocity  and  disturbance  velocity, 
respectively.  The  mean  flow  velocity  and  its  derivatives  are  computed  by  [1].  The  initial  acoustic  pulse  is 
given  by  (2.9).  The  uniform  grid  360  x  360  in  polar  computational  space  is  used. 

In  Fig.  13,  isobars  of  acoustic  pressure  at  i  =  0.5  are  presented  for  upstream  propagation  of  acoustic 
pulse  (Case  CYL-A)  at  M  =  0.5,  downstream  propagation  of  acoustic  pulse  (Case  CYL-B)  at  M  =  0.4,  and 
static  ambient  conditions  (Case  CYL-C).  In  Fig.  14,  the  centerline  pressure  distribution  is  presented  for 
these  cases  at  time  moments  t  =  0.5, 1.0, 1.5, 2.  In  comparison  with  the  Case  CYL-C,  the  acoustic  pressure 
is  amplified  at  the  neighborhood  of  the  centerline  in  the  Case  CYL-A  and  is  weakened  in  the  Case  CYL-B. 

To  get  the  sound  directivity  at  different  Mach  numbers  for  Cases  CYL-A  and  CYL-B,  the  root  mean 
square  of  acoustic  pressure  is  calculated  by 


Prms  — 


tfdt/T, 


(4.2) 


where  p  is  the  acoustic  pressure  (p  is  equal  to  p )  and  T  is  the  time  period  of  summation.  This  time 
period  is  chosen  constant  and  equal  for  all  numerical  experiments  presented  below.  For  temporal  numerical 
integration,  the  integral  in  (4.2)  is  computed  as  sum  of  values  of  squares  of  acoustic  pressure  at  all  time  steps 
over  the  period,  T.  The  p2  as  a  function  of  time  for  Cases  CYL-A,  CYL-B,  and  CYL-C  at  the  centerline 
location  -5Rcyi  is  presented  in  Fig.  15.  Curves  have  two  stairs  corresponding  to  passing  transmitted  and 
reflected  waves.  The  time  period,  T,  covers  the  pass  of  both  waves.  The  prms  represents  the  time-averaged 
potential  part  of  acoustic  energy  passing  through  a  fixed  point  in  space. 

In  Fig.  16,  the  prms  is  shown  as  function  of  angle  a  from  the  centerline  while  the  radial  coordinate  is 
equal  to  -5Rcyi.  The  curves  are  presented  for  various  Mach  numbers  from  M  =  -0.4  to  M  =  0.5  with  step 
0.1.  The  pressure  pulse  is  located  at  the  distance  0.25f?cy,  from  the  cylinder  (Fig.  16a, b)  and  at  the  distance 
0.5 RCyi  (Fig.  16c, d)  at  t  =  0.  The  positive  Mach  numbers  correspond  to  the  Case  CYL-A  (Fig.  16a, c),  the 
negative  Mach  numbers  represent  the  Case  CYL-B  (Fig.  16b, d),  and  the  zero  Mach  number  is  the  Case 
CYL-C.  The  latter  case  is  presented  in  (Fig.  16a-d)  as  a  benchmark  for  comparison.  The  presented  results 
are  normalized  in  such  a  way  that  the  pTms  at  the  centerline  in  the  Case  CYL-C  is  equal  to  unity. 

The  prTns  increases  with  the  Mach  number  and  reaches  2.2  while  M  =  0.5  (Fig.  16a).  Even  far  apart 
from  the  centerline  (a  =  60°),  the  prm»  for  M  =  0.5  is  about  twice  as  many  as  that  for  the  Case  CYL-C 
(M  —  0)  in  the  same  direction. 
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For  the  Case  CYL-B,  the  prm8  reduces  with  the  absolute  value  of  the  Mach  number  and  becomes  0.65 
at  the  centerline  while  M  =  —0.4.  The  level  of  prms  is  smaller  than  that  for  the  Case  CYL-C  up  to  a  =  45°. 

The  degree  of  increase  and  reduction  of  acoustic  pressure  is  slightly  bigger  when  the  pulse  is  located  more 
close  to  wall  (compare  Fig.  16a, b  with  Fig.  16b, d).  Note,  that  the  difference  between  acoustic  pressures 
is  quite  small  for  these  pulse  locations  in  Case  CYL-C.  The  reason  for  such  a  difference  for  non-zero  Mach 
numbers  is  that  the  reflected  wave  propagates  in  the  direction  opposite  to  the  far-field  sound  before  it  hits 
the  wall.  For  instance,  in  the  Case  CYL-A  the  reflected  wave  propagates  downstream  before  its  reflection 
from  the  cylinder  and  its  acoustic  pressure  decreases. 

5.  Conclusions.  The  behavior  of  acoustic  waves  originating  from  the  single  pulse  and  propagated  in 
the  non-uniform  incompressible  background  flow  is  studied  by  numerical  simulation  using  Hardin  and  Pope 
approach.  Higher-order  compact  spatial  finite  differences  and  Runge-Kutta  temporal  integration  are  used. 
Results  are  compared  to  those  obtained  for  the  same  geometry  in  the  static  ambient  conditions. 

It  is  shown,  that  the  acoustic  pressure  is  almost  doubled  while  the  waves  propagate  upstream  of  the 
stagnation  flow  and  the  background  flow  reaches  M  =  0.5.  Alternation  of  the  direction  of  stagnation  flow 
change  the  direction  of  sound  propagation  from  perpendicular  to  a  wall  to  parallel  to  it.  Sound  propagation 
in  a  corner  while  in  the  presence  of  the  background  inviscid  flow  leads  to  amplification  of  sound  and  flattening 
of  acoustic  waves. 

In  terms  of  acoustic  energy,  the  modified  sound  directivity  in  the  presence  of  stagnation  flow  is  mainly 
caused  by  redistribution  of  potential  and  kinetic  components  of  acoustic  energy  (while  sound  propagates 
upstream)  and  by  angular  redistribution  of  acoustic  energy  (while  directions  of  the  mean  flow  and  sound 
propagation  are  far  from  collinear).  The  pump  of  acoustic  energy  from  the  background  flow  by  means  of 
baroclinically  generated  vorticity  is  minor. 

The  propagation  of  acoustic  wave  originated  from  spherical  pulse  in  3-D  stagnation  flows  was  considered 
and  the  wave  pattern  in  the  neighborhood  of  the  single  stagnation  point  and  the  stagnation  line  are  discussed. 
The  amplification  of  sound  by  a  mean  flow  is  more  prominent  in  the  3-D  case  than  in  the  2-D  case  because 
of  steeper  acoustic  wave  bounds. 

Finally,  the  propagation  of  acoustic  pulse  in  the  flow  around  the  2-D  circular  cylinder  is  modeled. 
The  time-averaged  root  mean  square  of  acoustic  pressure  (rms)  is  presented  as  function  of  angle  from  the 
centerline.  To  compare  with  the  static  ambient  conditions,  the  rms  of  acoustic  pressure  at  the  centerline 
is  approximately  doubled  (upstream  pulse  position,  M  =  0.5)  and  halved  (downstream  pulse  position, 
M  —  0.4).  The  amplification  (upstream  pulse  position)  and  weakening  (downstream  pulse  position)  of  sound 
holds  for  a  large  angular  sector  apart  from  the  centerline. 

This  study  can  be  expanded  to  cases  where  a  compressible  mean  flow  is  computed  by  a  CFD  method  and 
the  acoustic  source  distribution  is  extracted  from  appropriate  turbulence  modeling.  Investigation  of  sound 
amplification  is  expected  to  be  important  for  higher  Mach  number  flows  where  the  velocity  and  pressure 
gradients  in  the  mean  flow  are  larger  than  those  in  the  low-Mach-number  flows  considered  here. 
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Appendix  A.  Mean-flow  variables  for  a  flow  between  two  straight  streamlines. 

In  Cartesian  coordinates,  the  stream-function  is  given  by 


=  A(x2  +  y2)05nsin(n  arctan(ar/2/)),  (A.l) 

The  velocity  field  (U(x,y)iV(x,y))  is  computed  by  differentiation  of  the  stream-function  defined  above 

(A.2) 

Using  Mathematica  [18],  the  following  expressions  for  velocities  and  their  derivatives  are  obtained 


U{x,y)  —  Anrn  2  (x  cos  nO  +  y  sin  nO) , 

V ( x ,  y)  =  Anrn~2  (y  cos  n6  -  x  sin  nO ) , 

—  Anrn~i((n  -  1  )(y2  -  x2)cosn9  +  2(n  -  I)z2/sinn0), 
ox 

v)  _  Anrn-4((n  _  -  y2) sinn0  -  2(n  —  l)xy cosnO), 

dy 


(A-3) 

(A.4) 

(A.5) 

(A.6) 


where  r  =  y/{x2  +  y2),  and  6  =  arctan^/y). 

For  an  incompressible  and  irrotational  mean  flow, 


dUjx, y)  _  dVj^y)  dUfay)  _  dV(x, y) 
dx  dy  ’  dy  dx 


(A.7) 
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Fig.  5.  Generation  of  acoustic  energy  ( Case  2-A):  a)  generation  of  vorticity,  and  b)  generation  of  acoustic  energy. 


Fig.  6.  Acoustic  energy  and  its  source  along  the  bisector  a  =  45°  at  t  —  0.5,  where  1)  is  the  acoustic  energy  (solid  line), 
and  2)  is  the  source  of  energy  (dashed  line).  L( 45°)  is  the  distance  from  the  origin  along  the  bisector 
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Fig,  7.  Acoustic  pulse  propagation  in  a  90°  comer  at  t  =  0.5  :  a)  in  presence  of  the  inviscid  flow ,  and  b)  in  the  static 
ambient  conditions . 


Fig.  8.  3-D  acoustic  pulse  propagation  in  the  static  ambient  conditions  in  presence  of  rigid  wall  Computational  domain: 
[-1, 1]  x  [-1,1]  x  [0,1]  1)  80  x  80;  2)  120  x  120;  and  3)  analytical  solution . 
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Fig.  11.  Acoustic  pulse  propagation  and  reflection  in  static  ambient  conditions  (analytic  solution):  a)  centerline  acoustic 
pressure ,  b)  the  spatial  derivative  of  acoustic  pressure  along  the  centerline.  The  direction  of  the  x— coordinate  coincides  with 
the  centerline,  graphs  ( 1 )  and  (2)  denote  2-D  cylindrical  pulse  (Case  2-B)  and  3-D  spherical  pulse  (Case  3-C),  respectively. 


21 


Fig.  12.  Potential  non-lifting  flow  over  a  circular  cylinder.  Here  P  denote  initial  position  of  acoustic  pulse ,  and  thick 
lines  denote  boundaries  of  cylindrical  computational  domain .  Isolines  of  mean  pressure  and  streamlines  are  shown . 
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Fig.  13.  Isobars  of  acoustic  pressure  at  t  =  1.5.  a)  upstream  propagation  of  sound  (Case  CYL-A)  b)  downstream  prop¬ 
agation  of  sound  (Case  CYL-B),  and  c)  static  ambient  conditions  (Case  CYL-C).  The  direction  of  x-  coordinate  coincides 
with  the  centerline .  Time  moments:  1)  t  =  0.5;  2)  t  =  1.0;  3)  t  =  1.5;  4)  t  =  2-0. 
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b  d 

Fig.  16.  Root  mean  square  (rms)  of  acoustic  pressure:  a,b)  pulse  location  Rinit  =  0.25 Rcyi,  c,d)  pulse  location  Rinit  — 
0.5 Rcyi,  a,c)  upstream  propagation  (Case  CYL-A),  b,d)  downstream  propagation  (Case  CYL-B) 
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